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ABSTRACT 

Observations of star forming environments revealed that the abundances of 
some deuterated interstellar molecules are markedly larger than the cosmic D/H 

■^ I ratio of 10~^. Possible reasons for this pointed to grain surface chemistry. How- 

ever, organic molecules and water, which are both ice constituents, do not enjoy 
the same deuteration. For example, deuterated formaldehyde is very abundant 

C^^ I in comets and star forming regions, while deuterated water rarely is. In this 

article, we explain this selective deuteration by following the formation of ices 
(using the rate equation method) in translucent clouds, as well as their evolu- 
tion as the cloud collapses to form a star. Ices start with the deposition of gas 

/\ ' phase CO and O onto dust grains. While reaction of oxygen with atoms (H or 

c3 I D) or molecules (H2) yields H2O (HDO), CO only reacts with atoms (H and D) 

to form II2CO (HDCO, D2CO). As a result, the deuteration of formaldehyde 
is sensitive to the gas D/H ratio as the cloud undergoes gravitational collapse, 
while the deuteration of water strongly depends on the dust temperature at the 
time of ice formation. These results reproduce well the deuterium fractionation 
of formaldehyde observed in comets and star forming regions and can explain the 
wide spread of deuterium fractionation of water observed in these environments. 

Subject headings: dust, extinction — ISM: abundances — ISM: molecules — 
stars: formation 
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Introduction 



Stars form through the collapse of interstellar clouds which are composed of gas and 
dust. As the cloud collapses, its density increases, its temperature decreases and it becomes 
shielded from external UV radiation that would otherwise dissociate chemical species present 
in the medium. In these regions, dust grains grow thick icy mantles from the deposition of 
atomic and molecular species onto their surfaces. This has been confirmed by observations of 
starless cores (Tafalla et al. 2006), which show that most gas phase species suffer a significant 
drop towards the core center. The missing gas species constitute the icy mantles that cover 
dust grains. As a star forms and heats up its environment, a portion of the icy mantles are 
released into the gas phase. This phase of star formation, called the hot core (massive stars) 
or hot corino (low mass stars) phase, exhibits a very complex chemistry rich in oxygen and 
nitrogen bearing molecules (Cazaux et al. 2003), but also shows species such as H2CO and 
methanol, which are building blocks for more complex organic molecules essential for the 
formation of life. Also, while deuterium is known to be 10^ times less abundant than hydrogen 
(Linsky 2003), many of these environments present a chemistry rich in highly deuterated 
species. In particular, deuterated methanol has been found in molecular clouds (Parise et al. 
2004) and in low mass protostars, with a maximum of ~ 30% for CH2DOH/CH3OH (Parise 
et al. 2002) and of ~ 3% for CD3OH/ CH3OH (Parise et al. 2004), while no detection has 
been reported in comets (Crovisier et al. 2004). Formaldehyde present a high deuterium 
fractionation as well (see Table 1), HDCO/H2CO is of a few percent towards protostellar 
cores (Roberts & Millar 2007, Bacmann et al. 2002) and towards hot cores (Turner 1990) 
and of a tenth of a percent towards hot corinos (Parise et al. 2006). Its doubly deuterated 
form D2CO/H2CO can reach up to 5% in hot corinos (Parise et al. 2006). Species that have 
been observed with 3 deuterium atoms (ND3; van der Tak et al. 2002) and CD3OII (Parise 
et al. 2004) are expected to be 10^^ times less abundant than their hydrogen form, with 
respect to elemental abundances. Yet, their enhancement can be of 12 orders of magnitude 
(Ceccarelli et al. 2007). 

This high degree of deuteration is attributed to grain surface chemistry. However, while 
water co-exists with methanol and formaldehyde in ices, it rarely has the same deuteration. 
In general, the degree of deuteration of organic material is very high compared to that 
measured in water (Ehrenfreund et al. 2002). The deuterium fractionation of water (Table 
1) can be about HDO/H2O~0.01% towards massive hot cores (Gensheimer et al. 1996) and 
~0. 02-0. 05% in comets and asteroids (Altwegg & Bockelee-Morvan 2003). These values are 
similar to SMOW (Standard Mean Ocean Water), which reports a fraction IIDO/H2O of 
0.015% in our oceans. However, some detections showed a high HDO/H2O ratio of few % 
towards the hot corino NGC1333-IRAS2A (Liu et al. 2011), and IRAS16283-2422 (Coutens 
et al. in prep., Parise et al. 2005). In this study, we address the selective deuteration of 
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water compared to formaldehyde. We used H2CO as a representative of organic molecules 
(H2CO is a building block of most of the complex organic molecules found in space), but 
our study also apply to other species formed with CO, such as CH3OH. First, we determine 
the chemical gas phase content of the region where ices form. Then, we use a rate equation 
method that follows the accretion and reactions of gas phase species that yield ices. Finally, 
we simulate how the ice compositions change as they undergo gravitational collapse to form 
stars. These latest results are compared with observations of star forming environments as 
well as cometary material. 



2. Origin of the ices 

Ices are formed in translucent clouds (clouds where H2 dominates over HI and carbon 
is converted from C II to C I to CO), in regions that are shielded from UV radiation. At 
extinctions Av>3 mag, species from the gas phase accrete onto dust grains and initiate the 
formation of ices (Whittet et al. 2001, Chiar et al. 2011). Therefore, the chemical content of 
the gas phase at this extinction is crucial to determine which species accrete onto dust and 
what the composition of ices is. Gas phase chemistry of dense cores (Roberts et al. 2003) 
shows that a high D/H ratio can be reached because of the inclusion of ions such as HD^ 
and D^ in gas phase chemistry, which produce deuterons upon dissociative recombination. 
This gas phase model is able to reproduce a high deuteration of formaldehyde in dense cores. 
Also, the degree of deuteration of molecular ions are sensitive to the orthoipara H2 ratio and 
hence to the chemical and thermal history of the precursor molecular cloud (Flower et al. 
2006). 

In this study, we concentrate on the ice formation driven by the accretion of gas phase 
species. To this effect, we determine the gas phase composition (in hydrogen and deuterium) 
of a medium submitted to UV radiation from neighboring stars, by using a model that 
couples gas phase and grain surface chemistry (Cazaux & Spaans 2009). This model follows 
the gas phase densities of molecular and atomic hydrogen as well as deuterium, as the 
medium becomes shielded to the external UV field (Draine & Bertoldi 1996). The results 
reported in figure 1 show the chemical composition in hydrogen and deuterium as function 
of the extinction, and therefore provide the initial conditions that are used in the next 
section (table 2) for the formation of ices (at Av~3). Note that neither oxygen, CO and ice 
chemistry are treated here, our goal being to assess the hydrogen and deuterium content of a 
medium in the region where ices are forming. As the hydrogen becomes molecular, deuterium 
can remain atomic, which leads to an enhanced D/H ratio at low extinctions. Then, when 
the medium is shielded, the chemistry lead by ions drives a D/H ratio that can reach 0.3 
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(Roberts et al. 2003). In Figure 1, the H/H2 and D/HD fronts are presented, as well as the 
D/H ratio as function of the extinction for an UV radiation field similar to the one in our 
local interstellar medium (defined as Go=l; Habing 1968). The D/H ratio increases with 
density for extinctions higher than 1 mag. In this study we consider the formation of ices in 
translucent clouds with initial gas phase conditions listed in Table 2. 



3. Dust favors deuteration 

In our grain surface chemistry model, we consider the formation of ices by accreting 
species from the gas phase, and let the icy mantles grow until the CO and oxygen from the 
gas phase are strongly depleted onto dust. Jeans instability, turbulence and/or triggering due 
to winds from nearby young stars will then lead to local density increase and gravitational 
collapse of a cloud fragment. During the collapse, the molecules remaining in the gas can 
still interfere with the ices and change their composition. Our model is therefore separated 
into two parts: 1) the slow formation of ices in diffuse environments; 2) the gravitational 
collapse of dense icy regions of the cloud beginning star formation. 



3.1. Translucent clouds 

3.1.1. Gas phase composition 

Ices form in shielded regions (Ay~3 mag), with CO and oxygen equally abundant in the 
gas phase (Hollenbach et al. 2009) and scale as 1.5x 10^^ n//, where n^^ the total hydrogen 
density. We assume that ices are first formed in translucent clouds, with initial gas phase 
conditions listed in Table 2. As time passes, oxygen and CO will mostly be depleted onto 
dust, impoverishing the gas phase. Their densities can be calculated as: 

—— = 1.5 X IQ-'^Uh - Odust X ^UdustCmsite, (1) 

1.5 X IQ-^^uh - C Odust X '^ndustcmsite, (2) 



dnco _T ,,,. Tn-4 



dt 

where Odust and COdust are the amount of oxygen and CO on the dust, in monolayers, 
nsites ~10^^ cm~^ is the density of sites, Tidust the density of dust (cm"^) and a their cross 
section (cm~^). Here we consider a grain size distribution that takes into account very small 
grain particles (Weingartner & Draine 2001), and derive ndMstO"/n_H'=2.8 10^^^ cm^. With 
this, we compute the densities of gas phase CO and oxygen with time, and also follow in 
parallel the chemistry occurring on the dust surfaces. 
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3.1.2. Grain surface chemistry 

The gas phase composition of the environments where ices originate, discussed above, 
is reported in table 2. We follow the population of species on the dust by using the rate 
equation method. Species present in the gas accrete onto dust at a rate: 

Race, = monolayer s~\ (3) 

where Uj and Vj are, respectively, the densities and velocities of the species i, and S is 
the sticking coefficient of the species with the dust (here we consider S=l). Since our rate 
equation model follows the populations of the species of the surface (with 1 monolayer =100% 
coverage), this rate is expressed in monolayer s~^. The binding energies between the chemical 
species and the grain surface are assumed to be weak (Van der Waals or physisorption), which 
is typical for icy surfaces. These energies are taken from previous studies (Cuppen & Herbst 
2007), with the exception of oxygen, that we take as 1390 K (Bergeron et al. 2008) and H2 as 
520K (Dulieu et al. 2005). The binding energies of CO and hydrogenated/deuterated forms 
are all considered to be 1200 K (Allouche et al. 1998). 

The species i present on the surface may return to the gas phase if they evaporate with 

Revap(i} =l^iX exp (---^), (4) 

Kb I 
where z/j is the oscillation factor of the species i (which is of 10^^ s~^ in physisorbed sites) and 
Ej is the binding energy of the species i, as reported in Cazaux et al. (2010) and discussed 
above. 

On the surface, atoms and molecules can travel on the dust with a mobility (Cazaux & 
Tielens 2004, Cazaux et al. 2010): 

E 

Ri = UiX exp{--^) (5) 

where Epp is the energy of the barrier between two physisorbed sites (considered as 2/3 of 
the binding energies Ej). 

Species on the surface can be photodissociated by UV photons (Cuppen & Herbst 2007), 
and can meet another species to form a product that will either stay on the surface or be 
released on the gas phase (Cazaux et al. 2010). In this study, we follow the coverage of the 
different species in order to determine the ices composition. 

We have constructed a chemical network that considers the hydrogenation and deutera- 
tion, but also the H-D and D-H exchanges for H2CO, HDCO and D2CO. The most important 
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reactions of the chemical network involving O and CO are presented in Figure 2 with the 
associated barriers for the reactions to occur, E^. The characteristic time for a reaction can 
be calculated as ~ 10^^ x exp(— 0.406 x ^/jired x -Eb) (oscillation factor times the probability 
of tunneling through a barrier of lA) , with fired the reduced mass for the reaction. This 
reaction occurs if its characteristic time is smaller than the time for the species to leave the 
site. 

The chemistry that involves oxygen has been described in a previous work (Cazaux 
et al. 2010). For the chemistry of CO, the different barriers for the reactions considered 
in our model have been derived by several authors (Fuchs et al. 2009, Hidaka et al. 2007, 
Hidaka et al. 2009). The reactions H2 + CO -> H2CO as well as HCO + H2 ^ H2CO + 
H, have very high barriers of 38400 K and 9000 K (Manion et al. 2008), respectively. Also, 
the reactions involving D2 are negligible (Hidaka et al. 2009). Therefore, we do not include 
reactions involving H2, HD and D2 in the CO chemistry network. 

Recent studies (Ratajczak et al. 2009, 2011) examined exchanges between H and D in 
ices composed of water and deuterated methanol as these ices warm up. D/H exchanges 
concerns the hydroxyl functional group of methanol and therefore occur between CD3OD 
and H2O that lead to CD3OH and HDO. In the ices, the amount of CD3OD is very low 
compared to the amount of water. These H/D exchanges can be responsible for a very small 
enhancement of HDO and are therefore not considered in this study. 



3.2. Collapse and star formation 

In the second phase of our model, the dust grains covered by ices are present in a denser 
environments that undergoes gravitational collapse. We assume the density of the medium 
to be nj:/~10^ cm~^ and our model cloud to collapse at the free-fall rate. The simulations 
follow the gas density evolution as: 

dnn uh ,„. 

dt tff 

where tff=y^37r/32Gp ~ 10^ yr is the free-fall time, with G the gravitational constant and 
p the mass density. Because CO and oxygen are mostly depleted onto the dust, only H, D, 
H2 and HD can accrete onto the icy surfaces with an atomic D/H ratio ~ 5 10^^. The initial 
gas phase composition of this step of our simulation is reported in Table 2. As the cloud 
collapses, and the medium becomes denser, the D/H ratio increases as shown in Figure 1. 
Molecular hydrogen, which is by far the most abundant molecule in the medium, accretes 
onto the last layer of ice. Once one layer of H2 covers the ice (around n//~10^cm~^), we 
assume that no further accretion is possible since the binding energies of species with H2 
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are very low, and therefore species bounce back into the gas phase. The ice constituents are 
followed during the collapse of the cloud, as a function of the density of the medium. 



4. Results 

We have computed the abundance of water and formaldehyde as well as their deuterated 
forms in interstellar ices. 

The oxygen accreting onto the dust transforms into water ice. Fig 3a shows the forma- 
tion of H2O and HDO on the dust surface as a function of time. The deuteration of water is 
very sensitive to the dust temperature. At low Tdust (12 K), H2 and H are both present on 
the ices, and oxygen can react both with H and H2. Because H2 is much more abundant than 
H, water forms through the reaction O+H2 while the formation of HDO involves atomic D. 
Therefore, the ratio HDO/H2O (fig 3b) scales with D/H2. At Tfiust= 15 K, short time scales 
allow a important amount of H atoms to be ice constituents, and oxygen reacts with atomic 
H to form water. As time proceeds, H atoms disappear by reacting with species present 
in the ices and H2 becomes largely dominant. Oxygen atoms coming from the gas phase 
then react predominantly with H2 molecules. Therefore, the ratio HDO/H2O scales with 
D/H for short time scales and with D/H2 for long time scales. At Tdust=^7 K, H2 molecules 
evaporate too fast to allow the reaction O+H2 to occur. Water is formed through successive 
hydrogenation of oxygen, and the ratio HDO/H2O scales with D/H. 

The CO accreting onto the dust transforms into formaldehyde. Fig. 3c shows the build 
up of H2CO, HDCO and D2CO ices. The formation of formaldehyde is made through the 
reactions CO+H— 7>HC0 and HCO+H— 7>H2CO. The CO-ices build up slowly, and CO present 
on the dust associates with an incoming H atom, or with a D atom. Therefore, HDCO/H2CO 
and D2CO/H2CO (fig 3d) scale with D/H and (D/H)2, respectively. 

Once the medium becomes denser, the cloud made of gas and dust (covered by ices) can 
undergo gravitational collapse and form a star. This star heats up its environment and the 
content of the ices is released into the gas phase (hot core/corino). At later stages, the star 
is surrounded by a disk, and the content of the ices can be present in the warmer parts, or 
can be locked into cometary material. In these different phases of star formation, deuterated 
water and formaldehyde have been observed. 

Deuterium fractionation of formaldehyde and water has been observed toward several 
astrophysical objects, as reported in table 1. While HDCO/H2CO and D2CO/H2CO ratios 
range between 0.01-0.28 and 0.01-0.04, respectively, HDO/H2O ratio are much more spread, 
ranging from few 10~^ to 0.02. In Fig 3e and Fig 3f, we compare the results of our model to 



these observations. 

Fig. 3f, presents the evolution of the deuterated forms of formaldehyde as the cloud col- 
lapses. As the ices are present in a denser environment with a higher D/H ratio, the fraction 
of deuterated formaldehyde increases by 1 order of magnitude. Our model reproduces the 
high degree of deuteration of formaldehyde as the cloud collapses. Also, the fractionation of 
formaldehyde is not sensitive to the temperature at which the ices formed in the translucent 
clouds. 

For water, on the other hand. Fig. 3e shows that the ratio HDO/H2O is constant as the 
cloud collapses. Once H2O and HDO are formed, a very important barrier (9600K) has to 
be overcome to break these molecules. Therefore, the deuterium fractionation is set during 
the formation of ices, in the translucent clouds, and deuterated water retains the memory 
of the ice formation. The ratio HDO/H2O is extremely sensitive to the dust temperature 
at which the ices formed in the translucent clouds. At low temperatures, water is formed 
with 0+ H2, while at higher temperatures H2 evaporate, and water is mode through O+H. 
As a result, our model can explain the spread in the HDO/H2O ratio observed in different 
astrophysical objects. 

The mantles species that are released in the gas phase could be transformed in other 
species through gas phase reactions. Goumans (2011) showed that gas phase reactions of 
H2CO with H(D) would lead to HCO + H2(HD), while H +D2CO would lead to CHD2O. In 
this case, the deuterium fractionation of formaldehyde could be increased by further reactions 
in the gas phase. However, Charnley (1997) show that the HDO/H2O ratio stays constant 
after ices evaporate, while HDCO/H2CO could be increased by a factor 5. Therefore, the 
deuterium fractionation found in this study, that is the result of the chemistry occurring 
in the ices, should not be dramatically changed by the gas phase chemistry that follows ice 
evaporation. 



5. Conclusions 

We used the rate equation method to simulate the formation of water and formalde- 
hyde, as well as their isotopologues on dust grain surfaces. We model the formation of 
ices in translucent clouds and incorporate these ices in denser environments that undergo 
gravitational collapse to form stars. 

We find that HDCO/H2CO and D2CO/H2CO depend on the gas D/H ratio which 
increases as the medium collapses and becomes denser. 
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For water, on the other hand, HDO/H2O strongly depends on the dust temperatures at 
which ices form since water formation involves H2 at low T^ust (H2 very abundant on dust 
surface) and H at higher T^ust (H2 evaporate). Consequently, the degree of water deuteration 
is directly related to the physical conditions (i.e., dust temperature) at which ices form in the 
quiescent cloud just before it collapses to form stars. H2CO is the simplest representative of 
organic species formed via hydrogenation of CO onto dust surfaces. Our results can then be 
generalized to other species (such as CH3CO), which also depends on surface hydrogenation 
of CO. This explains the selective deuteration of water and organic molecules observed in 
the ISM, similar to that found in comets, IDPs and chondrites (Ehrenfreund et al. 2002, Fig. 
9). 
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Table 1: Deuterium fractionation of formaldehyde and water observed in different environ- 
ments 



Environment 
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'^Hershel results IRAS16293-2422 Coutens et al. in prep. 
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Fig. 1. — Left panel: atomic (H and D) and molecular (H2 and HD) hydrogen and deuterium 
abundances with respect to total hydrogen (nHI+nH2) in a cloud with a density of 10 
(dashed), 100 (solid) and 10^ cm^^ (dotted) as a function of the visual extinction (mag). 
Right panel: Variation of the D/H ratio as a function of the visual extinction (mag) for 
densities varying from 10 to 10^ cm^^. The line shows Ay~3, where the ices initiate. 
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Table 2: Initial conditions of translucent (at Ay~3) and collapsing clouds used in our simu- 
lations 
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Fig. 2. — Principal reactions involving CO (left) and O (right) that are considered in our 
network. The reactions in boldface are associated with barriers (in Kelvin). Barriers that 
have not been determined, but that we consider as the mean of similar reaction are denoted 
with an asterisk (*). Note that O reacts with atoms and molecules, while CO only reacts 
with atoms. 
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Fig. 3.— H2O, HDO, H2CO, HDCO and D2CO content on the dust during (1) the formation 
of ices in a translucent cloud and (2) the gravitational collapse of a cloud to form a star. 
Lines represent "Y^dust =17 K (solid), 15 K (dotted) and 12 K (dashed). Top: surface coverage 
in monolayers (%) as function of time for H2O and HDO (a) H2CO, HDCO and D2CO 
(b). Middle: deuterium fractionation of water (c) and formaldehyde (d). Bottom panels: 
model versus observations of HDO/H2O (e) HDCO/H2CO (filled symbols) D2CO/H2CO 
(open symbols) (f) towards several astrophysical objects as a function of total hydrogen 
density (Table 1). 



